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Abstract 

In this paper, we show that the kernel function of Cauchy type for type BC intertwines the 
commuting family of van Diejen's (/-difference operators. This result gives rise to a transforma- 
tion formula for certain multiple basic hypergeometric series of type BC. We also construct a 
new infinite family of commuting q-difference operators for which the Koornwinder polynomials 
are joint eigenfunctions. 
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1 Introduction 

In the theory of Macdonald polynomials of type A, the kernel function of Cauchy type has been 
used to derive various important properties of Macdonald polynomials [TT] I14j . Kajihara's 
Euler transformation formula for multiple basic hypergeometric series can also be regarded as an 
application of the kernel function of Cauchy type [6]. 
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Recently, Y. Komori, M. Noumi and J. Shiraishi in [8J introduced the kernel function <£(x; y\q, t) 
of type BC in the variables x = (xi, . . . , x m ) and y = (yi, . . . , y n ) relevant to Koornwinder poly- 
nomials. The kernel function <&(x;y\q,t) satisfies the following g-difference equation: 

(t)DfO>(x;y\q,t) - (t)D^(x;y\q,t) = {t m ){r n ){a 2 t m ~ n - ^(x; y\q, t), (1.1) 



where a = \Jabcdjq and (z) is the multiplicative notation for trigonometric function 

(z) = z^ — z~2 = — z~*(l — z). (1.2) 
In this identity, Df is the Koornwinder ^-difference operator in the x variables 

m m 

Df = Y J Mx)(T q , Xl - l) + Y,Mx- 1 ){T-X - 1), (1.3) 

i=l i=l 

{aXi}{bXi){cXi)(dXi) -r-r (tXjXj) (tXjXj 1 ) 

i[x) <*?>«> (^x^xj 1 ) ' 1 • j 

Tq^xif (%1 j • • • j Xj, . . . , X m ) — f(%l j • • • j (?Xj , • • • , X m ) , (l - ^) 

and denotes the Koornwinder operator in the ?/ variables with the parameters (a, 6, c, d) replaced 
by (y/tq/a,y/tq/b,y/tq/c,y/tq/d). In this paper, we show that $(x;y\q,t) intertwines the whole 
commuting family of van Diejen's g-difference operators, which includes the Koornwinder operator 
as the first member. 

In Subsection 12. 1\ we recall some basic facts on van Diejen's g-difference operators. We also 
state our main result in Subsection 12.21 and prove it in Subsection 12.31 In the proof of the main 
result, we show a rational function identity of x variables and y variables. By using a method 
of principal specialization, from this identity we derive two types of transformation formulas for 
multiple ^-series (Theorem 13.11 and Theorem 13. 3p in Section [3l Theorem 13.31 recovers one of the 
C type transformation formulas, due to H. Rosengren |19j . In Section U from the special case of 
Theorem 13.11 we construct a family of explicit q-difference operators of "row type" for which the 
Koornwinder polynomials are the eigenfunctions. 

Throughout the present paper, we assume that the base q is a complex number such that 
< \q\ < 1. We also assume that a, 6, c, d, q, t are generic complex numbers. We use the notation of 
fractional powers a K b x ■ ■ ■ x^ ■ ■ ■ of multiplicative variables a,b, . . . x, . . . as an conventional notation 
ex.p(nA + XB + • • • fix + • • • ) under a parametrization a = exp(^4), b = exp(_B), ... ,x = exp(X), . . . 
by additive variables A, B, . . . , X, 

2 Kernel identity of Cauchy type 
2.1 Van Diejen's g-difference operators 

In this subsection, we recall some basic properties of the family of van Diejen's g-difference opera- 
tors. For further details, we refer the reader to [2j HI [8j |9] . 
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The family of van Diejen's (/-difference operators {D^(a, b, c, d\q, i)}^=o i n the variables x is 
defined as follows: 

D% := D*(a,b,c,d\q,t) = £ V^^U^^x)^, (2.1) 

/C{l,...,m} 

0<|/|<r 
ei=±l(ie/) 

_ i r (<%&xf,cxf,(fa^) t r {tx^xj ,tqxfxj) (tx?xf) 

Uj, r (x)- 1^ ("1)11 (x 2S t -2SK 11 — OTT J-^- 11 T^ilT' (2 " 3) 

|/|=r i<j ieJ\/ 

5i=±i(ie/) 

where I c = {1, . . . , m} \ I and T q ^ = Yl ieI ^q%i- We also used the shorthand notation 

(zi, • • • ,z k ) = (zi) ■ ■ ■ (z k ), (zw ±:L } = {zw,zw~ l ). (2.4) 

We will use the following notation of (/-shifted factorial in this paper: 

l 

(*) q ,i = = (-l) l q-*®z-* (z; q) t (I = 0,1,2,...), (2.5) 

where (z;q)i = n!=i(l — q l ~ x z). For these two types of (/-shifted factorials, we use the shorthand 
notation as 



{zi,...,z k ) q ,i= H {zi) q ,i, {zw ±l ) q + = {zw) q ,i{zw l ) qh (2.6) 

l<i<fe 

(z 1 ,...,z k ;q)i= Yl (zf,q)l, (zw ±:L ;q)i = (zw;q)i(zw~ 1 ;q)i. (2.7) 



i<i<fc 

Let w(z) and v(z) denote the following rational functions, respectively: 

, . (az, bz, cz, dz) , . (tz) 

( ») ■ «w-y- < m > 

Then V € i : j(x), Uj^(x) are also expressed as 

V e i,j(x) = Y[w(x?) Jl v{xfx^)v(qxfxj) ]J v{x^xf), (2.9) 

iei i,jei iei 

i<j i&J 

uja*) = E c-irlH*?) II ^N^^V^') n t,(xf*xf ), (2.10) 

|/|=r i<j jeJ\/ 

Si=±i(iei) 

where v (xixf 1 ) means v(xiXj) ■ v{xiX^ j 1 ). 
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Let W m be the Weyl group of type BC m acting on the Laurent polynomials in the variables 
x = (x\ , . . . , x m ) through the permutations of the indices and the inversions of the variables. 
Under the assumption that a, b, c, d, q, t are generic, for each partition A = (Ai, . . . , A m ) there exists 
a unique Wm-invariant Laurent polynomial P X {x) = P\(x;a,b,c,d\q,t), called the Koornwinder 
polynomial attached to A, satisfying the following conditions [3]. 

(1) P\(x) is expanded by the orbit sums m Ai (x) = YlueW.a xV as 

P x (x) = m x (x) + ^2 c \^ m v( x ), ( 2 - n ) 

where c Xfl £ C and < means the dominance ordering of the partitions. 

(2) P\(x) is a joint eigenfunction of van Diejen's ^-difference operators D x : 

D x r P x {x) = P x (x)e r (atP™q x ; a\t), (2.12) 

where p m = (m— 1, . . . , 1, 0) and e r (x; a\t) are the interpolation polynomials of column type defined 
by 



e r (x;a\t)= J2 ; f 1_1 a)e(a; i2 ; f 2_2 a) • • • e(x ir ; f r " r a) (2.13) 

l<il <---<i r <m 

^2 e(x il ;t m - h - r+1 a)e(x i2] t m - i2 - r+2 a)---e(x ir -t m - i ^a), (2.14) 

l<il <---<i r <m 

e(z; w) = (zw)(zw~ 1 ) = z + z^ 1 — w — w~ . (2-15) 

Note that e r {x; a\t) is VFm-invariant and satisfies the following interpolation property (See [8]): For 
any partition /i 7$ (l r ), 

e r (at Pm q^;a\t) = 0. (2.16) 

2.2 Main result 

We recall the definition of the kernel function <3?(x; y\q, t) of Cauchy type associated with the root 
systems of type BC in the variables x = (x±, . . . ,x m ) and y = (y±, . . . ,y n ). The kernel function 
$(x;y\q,t) is defined as a solution of the following linear (/-difference equations: 

T q , Xi $(x;y\q,t) = $(x;y\q,t) TT e (v^M^) (l<i<m), (2.17) 

T qm <Z>(x;y\q,t) = <Z>(x;y\q,t) J] <^^A (\<k<n). (2.18) 
Such a <&(x\y\q,t) is a multiple of the function 

^{x;y\q,t) = {x x ...x m r \{ ^^ q) ~ , (2.19) 

l<k<n 

by a g-periodic function with respect to all the variables x and y. Here (z; q)oo = Il£o(^ ~~ 1 lz ) an( ^ 
7 is a complex number such that q 1 = t. We note that four types of explicit formulas for kernel 
function of Cauchy type including <&q(x; y\q, t) are introduced in [8]. 
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For any integer I, let e{z; w)gi be the q-shifted factorial of type BC with base point w defined 

by 

{e(z; w)e(z; qw) ■ ■ ■ e(z; q l ~ 1 w) (Z > 0), 
1 n <Q) (2-20) 
e(z; q l w)e(z; q l+1 w) ■ ■ ■ e(z; q^w) 

We also define a generating function of D x and that of D y r by 



V x {u) := V x (u- a, b, c, d\q, t) = £(-l) r L>*e(u; a) t , m _ r , (2.21) 

r=0 

n 

T>y{u) := Vy(u;a,b,c,d\q,t) = £(-l) r 5»e(«; 5) t , n _ r , (2.22) 

where (a,b,c,d) = (y/tq/a,y/tq/b, y/tq/c, y/tq/d), so 5 = i/a. We also denoted Z)^ (a, b, c, d\q, t) 
by Dr- For any function /(z) = f(z;a,b,c,d) depending on the parameters (a,b,c,d), we write 
f(z) = f(z;a,b,c,d). Then we have the following theorem. 

Theorem. 2.1. The kernel function $>(x;y\q,t) intertwines the q-difference operator D x (u) in the 
x variables with the q-difference operator T> v (u) in the y variables: 

V x (u)<f>(x; y\q, t) = e(u; a) t , m - n T> y {u)®{x- y\q, t). (2.23) 

We call this equation a kernel identity of Cauchy type. 

A proof of this theorem will be given in the next subsection. We now give some remarks related 
to Theorem 12.11 Firstly, it is known that Theorem 12.11 in the case of n = holds. Namely, the 
constant function 1 is the eigenfunction of van Diejen's ^-difference operators [3]: 

V x {u)-l = e{u-a)t,m. (2.24) 
We will use this fact as the starting point of our proof. It is also known by [8] that 

m m 

£(-l) r e r (x; a\t)e(u; a)t, m -r = e(u; Xi). (2.25) 

r=0 i=l 

In general, for any partition A we have 

m 

V x (u)P x (x)=P x {x)\{e{u-at m -*q^). (2.26) 
i=l 

Secondly, comparing the coefficient of e{u\ a)t )m -i in the left-hand side of f|2.23|) with that in 
the right-hand side, we obtain (11, lj) . In fact, the g-Saalschiitz sum gives the transformation formula 
for the base points of the ^-shifted factorials of type BC: 



e{w;b)t,i= £ (-1)' 



0<r<l 



(t§(l-i)6;t3(i-0 /a)t, r e(w;a)u- r , (2.27) 
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(-1)' 



<*-'>*.r 



(t) 



t.r 



(2.28) 



It follows from this formula that 

e(u; a)t ym -ne(u; t/a) t n - k 



E n ^ e(t^ n - k+1 ^/a-,th(^~^+ k )/a) t!l e(u-a) t ^ k ^ (2.29) 



0<l<n-k 

Comparing the coefficients of e{u\ a)t )m -r in the both sides of (|2.23p . we have 



D?$(x;y\q,t)=J2D V k 



k=0 



n — k 
r — k 



,^ 2 (n-k+i) /a . t !(l+»-2"*fc)/ a ) t r _ fc $(a. ; y \g } t). (2.30) 



The formula (|2.30|) for r = 1 recovers the result (jl.ip of [8]. 
2.3 Proof of the main result 

It is enough to show the case where m > n > 0. The identity (|2,23p is equivalent to 

<&{x; y\q, t)- l V x (u)<$>{x; y\q, t) = i)" 1 ^^^; t)e(u; a\ m -n- (2.31) 

Regarding this as a rational function identity of the variable y n , we prove it by computing the 
residues and the limits as y n — > oo. 

The generating function D x (u) is expanded as 



(2.32) 



/C{l,...,m} 
e i= ±l(ieJ) 



Uj(u;x)= ^2 e ^ a )t,\J\-\i\^w{x^)^v{xfx S P)v{q l x i Si x j 53 



icJ 
8i=±i(iei) 



i<j 



) n vi^xf). (2.33) 
j'eJ\/ 



Similarly, we expand D y (u). We also define the rational function F(z;w) in the variables z 
(zi, ... ,z r ) and w = (wi, . . . ,w s ) by 



F(z;w)= H 



l<i<r 
Kk<s 



e{\jq/tzj;w k ) 
e(^/tqzi;w k ) 



(2.34) 



For any subset I = {ii, . . . ,i r } C {1, . . . , m}, \I\ = r and signs e, = ±l(i G I), we write x\ 
(Xj* 1 , . . . , x-* r ). Then, (|2.3ip is expressed as 



2 (-l) |/| K/,/ C (x)C// c (u;x)F(x};y) 



IC{l,...,m} 

ei=±i(ie/) 



e(u;a) t 

,m— n 

E (_l)l^ly e ^ Xc (y)C/^ c ( n;y )i?(^ ;rE ). 



(2.35) 



K{l,..,n} 
e fc =±l(feeJC) 
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We prove this identity by induction on re, starting with (|2.24p of the case re = 0. We assume that 
our identity holds when the number of y variables is less than re. 

Firstly, we consider the residues of the both sides as meromorphic functions in y n . The left-hand 
side of (]2.35p may have the poles at 

y n = Vtqxf\ -^=xf l {\<i<m). (2.36) 

On the other hand, there may be the poles at 

y n =y/Tqxf\ -^xf 1 (1 < i < rre), ±1, ±q 1 /\ 

V*9 (2.37) 

±q~ 1/2 , vt\ QVk 1 , Q'^f 0-<k<n) 
in the right-hand side. However, we can check by direct calculation that the points other than 
Un = Vtqxf 1 1 —j=xf X (i = 1, . . . ,m) are apparent singular points. Since (|2.35p is invariant under 



ftq 

the inversions and permutations for x and y, we have only to analyze the residue at the point 

Dn = ^/tqXjn. 

In the left-hand side, the term indexed by (I, e) has a pole at y n = \/tqx m if and only if rre € / 
and e m = 1. Note that 

F(x],;VT q x m ) = [J ^ XmX ^ tXm/ , X i\ = II ( 2 -38) 

(tqXmX^Xm/x^) V^X^X^) 

where I' = I \ {m}. Thus it follows that the residue is equal to 



Vt(ax m , bx m , cx m , dx m )(y/tq- y/q/i)x 

T~2 + 2~\ J 1 v \ x j X m ) 

[x m ,iqx m ) l<j<m-l 

■ F(x m ;y') x (l.h.s. for the case of (x';y f )), (2.39) 

where x' = (xi, . . . ,sc m _i) and y 1 = (j/i, . . . ,y n -i). 

In the right-hand side, the term indexed by (K, e) cannot have a pole at y n = ^Jtqx m unless 
re G K and e n = —1. The corresponding residue is equal to 

, bx m , cx m , 

/ 2 + 2 \ 11 v \ x j X m ) 

[x m ,iqx m ) l<j<m-l 

■ F(x m ;y') x (r.h.s. for the case of (x';y')). (2.40) 

Therefore it follows from the induction hypothesis that the residues of the both sides at the point 
Vn = Vtqx m are equal. 

Next, we calculate the limits of the both sides as y n — > oo. It is easy to check 

lim (l.h.s.) = (l.h.s. for the case of (x;y r )). (2-41) 

We consider the limit of the individual terms of the right-hand side in the following three cases: 
(i) n G K and e n = 1, (ii) n G K and e n = —1, (hi) n ^ K. 
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By direct calculation, we can check in the case (i) and (ii) respectively as follows: 

( ^ 

Y,(-±) lKl VeKMy)UKc(u;y)F(y'i < ;x 

i n£K 
\e n =l 

_ at n~ m -i ^ (-l) lK ' l V eK/ ^(y')U K/ c( U] y')F(y^;x), 



lim 



J 



(2.42) 



K'c{l,...,n-1} 
e k =±l(k&K') 



( 



lim 



-a 



Y, (-l) lKl V €K My)UKc(u;y)F(y* K ;x) 

\ n&K 

\e»=-l 

l t m-n+l J- (-l)^V eKl!K/ c(y')U Kl c(u-,y')F(yh-,x). 



(2.43) 



K'c{l,...,n-1} 

e k =±l(k&K') 



In the case (iii), we divide the sum 



U K c(u;y)= Y e ( u ^ a )t,\K-\-\L\]Jw(yl k ) Y[ viy^y^viq 1 y k 5k y l Sl ) ] J v{y S k k yf l 



LcK c 
S k =±l(k&L) 



k,l&L 
k<l 



l€K c \L 



into the three groups of terms as 

(a) n G L and 5 n = 1, (b) n G L and 5 n = — 1, (c) n ^ L. 
Combining the limits of these three cases, we obtain 



(2.44) 



lim J2(- 1 ) lKlv ^My) U K<u;y)F(y e K ;x) 



(u + u- 1 ) V (-i)^V £ ^-(y / )^(^;2/ , )i ? (y^;^ 



(2.45) 



A"c{l,...,n-1} 
e k =±l(k€K') 

From (H32|), fl!P3"D and (i2Ti5l . it follows that 

lim (r.h.s.) = (r.h.s. for the case of (x;y)), (2.46) 

and hence we complete the proof of Theorem 12.11 

Replacing (q, t) by (t,q), the formula f|2.35|) can be rewritten explicitly as follows. 

Theorem. 2.2. Given two sets of variables x = (xi, . . . ,x m ) and y = (yi, . . . ,y n ), the following 
identity holds: 

\ \ > 11 ^ x 2ei £j,2e^ 11 l„H^ e i +„ei„ £ i\ 11 



IC{l,...,m} 

ei=±i(iei) 



s 



\- ( t f-r: v tt (axf ,bxf ,cxf,dxf) , r (qx*xf,tx?x//q) 

Jci c \ ieJ \ « ' x « ' ijeJ \ x i x j x j / 

5i=±i(ieJ) i<i 

A„.±i\ 



jGl c \J l<k<n 

Kc{l,...,n} \ fcex ^fe '^fc / 

e fe =±i(fcejsr) 

t-t {gy e k k y e l l ,tqy e k k y e l l ) y-r (qy^yf 1 ) v / , /-, , 

fe<i *eA" c 5 fc =±i(fcei) 

tt {Vtqy 5 k k /a, \/tqy 5 k /b, Vtgy s k k /c, Vtqy s k k /d) yr (gy s k k y s l l ,ty s k k y s l l /q} 

keL (y 2 k k ^y 2 k h ) k ,iGL (y 5 k k y s i l ^yl k yi l ) 

k<l 

1&K C \L l<i<m 



3 Transformation formulas for multiple basic hyper geometric se- 
ries 

In the case of type A, the kernel function of Cauchy type intertwines the Macdonald g-difference 
operators [16]. This property gives the rational function identity which is similar to (|2.47p . Applying 
certain specializations to this identity, Kajihara [6] derived the Euler transformation formula for 
multiple basic hypergeometric series. In the same way, we propose two types of transformation 
formulas for multiple basic hypergeometric series. 



3.1 Type BC case 

In this subsection, we derive a transformation formula of type BC from Theorem 12.21 We take 
the multi-indices a = {ax, . . . , «m) £ N Af and (3 = {fix, . . . , /3jv) £ such that \a\ := Yli=i a i = 
m J \P\ = Ylk=i @ k = n - Here N is the set of non-negative integers. Then we consider the following 
specializations: 

x = p a (z; q) := (zx,qzx, g" 1 " 1 ^; z 2 , qz 2 , q a2 ~ 1 z 2 ) ■ ■■\z M , qzM, q aM ~ 1 z M ), (3.1) 
y =p /3 (w;q) := {wx, qwx, ■ ■ ■ , q Pl ~ l wx; w 2 , qw 2 , ■ ■ ■ , q l32 ~ 1 w 2 ; ■ ■ ■ ; w N , qw N , . . . , q^ N ~ 1 w N ). (3.2) 

These specializations are called multiple principal specializations. We apply these to (|2.47p . 
For each pair (I, e), we divide the subset I as I + U I~ by setting 

1+ = {i£ I\ €i = 1}, r = {ie I\ €i = -1}. (3.3) 
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Similarly, we divide the subset J by J + = {i £ J\8{ = 1} and J = {i £ J\5i = — 1} for each pair 
(J, 5). We also denote the complement {1, . . . , m}\(I U J) by C. Using these symbols, we rewrite 
the left-hand side of ff2TTT|) as 

El i -i i i 71 \ TT (axi,bxi, cxi, dxi) 
(-1)1 1+1 \e(u;y/q/ta) qt{cl [[ 

I+UI-UJ+UJ-UC \ i€l+Uj+ V ll 

={!,.. .,m} 

tt (a/xi,b/xi,c/xi,d/xi) yr {qxiXj) t .2 rr (qx^xj 1 )^ yr {qxjXj 1 ) t ,2 

i<i »<3 ie-f 

-|-r (ggjgf 1 ) rr (g^r 1 ^ 1 ) TT (qXiXjitXjXj/q) yr {qx^j\tx^^^ 

ieJ _ ucuj+ j'eJ _ ucuj+ i<j i<j 

{qxjx- 1 ,tx iX ^ 1 / q) (qxjxf 1 ) (qx^xf 1 ) {y/t/qxjyf 1 ) yr {^/tf^jg^\ 

jgj- j'GC jeC l<k<n l<k<n 

(3.4) 

We first consider the principal specialization x% = q l ~ 1 z (1 < i < m) of a single block. Noting 
that (qxi/xi+i) =0 (i = 1, . . . , m — 1), we find that non-zero terms arise from the divisions of the 
following form: 

r = {l,2,...,n}, J- = {i 1 + l,i 1 + 2,...,i 2 }, C = {i 2 + l,i 2 + 2,...,i 3 }, 

J + = {is + 1, i 3 + 2, . . . , I + = {i 4 + 1, 14 + 2, . . . , m}, < i\ < i 2 < h < U < m. 

Next, we consider the multiple principal specializations x = p a (z;q),y = pp{w;q). We replace 
the index set {1, . . . , m} by 

{l,...,m} = {(i,o)|l <i<M, < a < a { } (3.6) 

and write xu^ = q a Z{. For any two multi-indices fi, v € N A:f , if /^j < i/j (i = 1, . . . , M) then we write 
fi < v. From the same argument as above applied to each block, I~, J~, C, J + ,I + are parametrized 
by the four multi- indices fx~, v~, fi + 6 N M such that < < V < u + < ji + < a as follows: 

I~ = {(t,o)|l < i < M, < a < ^r}, 
J" = {(i,o)|l < i < M,fir < a < 

C = {(t,o)|l < i < M, vr <a< v+}, (3.7) 
J+ = {(i,a)|l <i<M, vf < a < //+}, 
J+ = {(i,a)|l <i< M, fif <a< an}. 

In the following, we omit the base (7 in the g-shifted factorials (z) q> k and e(z; w)q t k- 

With this parametrization (|3.6p of indices, the formula fj3.4j) specialized by x = p a (z;q),y = 
Pp(w\q) gives rise to 



n (azi,bzi,czi,dzi) ai -pr {q aj ZjZj,tq a i ZjZj) ai yr {^/tfqziW kl s/tqzij l q l3k w k ) t 



1 , s/ ( z i, tz i)ai i<i<j<M (ziZjitZiZj)^ (q^ y /t/qz i w k ,y/tqz i /w k ) ai 

Kk<N 
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_ ( i i , i ,,,, _. , . „ (z l /a,z i /b,z i /c,Zi/d) - 

(_i)M+KM* le(u;v^A«) K |_ M [] 



(azi,bzi,czi,dzi),,+ 

0<fi-<u-<u+<fi+<a \ l<i<M X l ' " l ' 

"I r (q 1 *' ZjZj/tq,tq^ +m j ZjZj/q) -p-r (g^ ~ M J Zj/zj,q^ ~^ Zj/zj) 
l<kj<M (ziZj/tq^ZiZj/q) l<i<j<M <*AW*i) 

(g"* - "j Zi/Zjrfi-^ z i /z j )(z i Zj/t,z i Zj/q) u - +u - (q»i+4 ZiZj / q ) 

l<i<j<M ( Z i/ Z 3, Z i/ Z j)( Z i z j> tz i Z j/<l)v++v+ l<i}<M ( Z i Z i/l) 



(g"" +v 7 Zi Zj/q, g M « Zi/tz^iZiZj/q}^-^ (tZiZj /q) ^+ +u + (frZj/tq)^ (ztZj/q)^ 
UJ<M {ziZj/q, Zi/tqfi Zj) {z l z j /t) ^ +v - (z iZj ) ^+ (q a J z iZj ) ^- (tq a i z { Zj) ^ 



,+ 



n 



{Zi/q"i zj, Zi/tq^ Zj) -{zi/q a ^ zj) +(qzi/tq^ Zj,Zi/q^ Zj ) -{zi/q^ Zj) -t 



Ki 



,j<M {qzi/tq v i Zj,qzi/zj) -(qzi/zj) +(qzi/tq^ Zj,qZi/zj) -(qzi/zj) + 



n{q l3k z i w k / y /tq, JqJtZi/w k ) - (g & ^/t/qziW k , y/tqzi/w k ) + \ 
7= 7= « • ( 3 - 8 ) 
i<i<M ( Z i w k/Vtq, ^q/tZi/q^ Wk )^{^t/qziW k ,y/tqzi/qPk Wk )^+ J 

l<k<N 

We used {a) k version of some formulas of Appendix I in |5j. We also used the formula 

n (qx i /xj) Xl _ -i-r (q X ^Xj/xj) 
(x l /q^X i )y 11 (Xi/Xi) ' K } 

l<i^j<m X 11 H ■ ?/Al l<i<j<m N ' 31 

due to Milne [131 Lemma 6.11]. 

We apply the same specializations to the right-hand side of (|2.47|) . If we denote (|3.8p by 
(z; w; a, b, c, d), Theorem 12.21 implies the following duality transformation formula: 



.y J (z;w;o,6,c,d) = e(u; ^fq/ta) M -NF^ a \w]z;^/a,y/tq/b,^/c,y/tq/d). (3.10) 

Relabeling M and iV by m and n, and replacing the variables and parameters by 

Zi^y/tqxi(i = l,...,m), w k ->• y k (k = 1, . . . ,n), 
(a,6,c, d) -> {y/tq/a 1 ,^/tq/a 2 ,\/tq/az,\/tq/ai), 

we obtain the following theorem. 



Theorem. 3.1. Xei ao = y 01020304/5. Tafce £wo sefe 0/ variables x = (xi, . . . , x m ),y = 
(yi,...,y n ) and two multi-indices a = (ai,...,a m ) S N m ,/3 = (/3i, . . . , /3 n ) G N n . TTien i/ie 
following identity holds: 



l r {tqxj/ai , tqxi/a 2 , tqxi/a 3 , tqxj/a^ -p-r {tq a i +1 XiXj,t 2 q a i +1 XiX j) ai 
(tqxl^qxf)^ i<i<j<m {tqXiXj,t 2 qXiXj) at 



l<j<m xy iy,c ' y ' ,M,fl ' 0<u-<r<i/+<u+<a\ 
KJfc<n 
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iXj I 



{a 1 x i ,a 2 x i ,a 3 x i ,aiX i ) u - ^ +^ XiXjj t 2 q^ + ^ x, 

l<i.<m { t< l x il a l^ < l x il a 2,tqXi/a^,tqXi/a i ) v + 1<;i 4j <m (x i X j ,t 2 XiX j ) 

(q^ Xi/xj,q^ Xi/xj,q u * ~ v j Xi/xj,q"> ~ v i Xi/xj)(tXiXj,qxiXj) - , - 
J j 1— J- 

{x{/ Xj , Xij Xj , X{/xj , x%j Xj) (t^XiXj , tqx{Xj) u + \ u + 

i\?<,j\?n i 3 

n(tq^ XiXj,tq u * +u 3 XiXj,q^ Xi/txj)(txiXj) - , + {t 2 XiXj) + , + 

l<i, 7<m (tX{Xj, tX{Xj, Xi/tq^i x j}(q x i x j)a^+u^ j) u + +u~ 

n(xiXj,Xi/q^ Xj,Xi/tq a JXj) -(txiXj,Xi/q a ixj) +(qxi/tq"j Xj,Xi/q u i xj) - 
i<i,j<m (tq a 3 +l XiXj,qXi/tq v i Xj,qxi/xj) - {t 2 q a i +1 XiXj,qxi/xj) ^+ (qxi/tq^ x j ,qx i /x j ) l 
(xi/q^i Xj) v + 1 r (q^XiVk^qXi/yk) ^-(tqP k Xiy kl tqxi/y k ) ^+ \ 



rr lii i^L rr 

l<ti<m ^ Xi / X o)vf l<i< m ( Xiyk ' ql *Wl/*>„- ttf 1 Pkx i/yk)^J 
l<k<n 

= e(u- y/tq/ao) TT ( aiyk,a2yk ' a3yk,a4yk ^k TT kq^VkVlM^VkVl} fa 

n ^S V A e ((-D^'^'^-'e^^),.,-,-, 

i< fc <„ i v^yk/^)h <A-< K -< K +<A+<A 

l<i<m 

„ (yk/a 1 ,y k /a 2 ,y k /a 3 ,y k /a A ) K - ( q ^+K y kyi /tq,tq x t +x t y k y x /q) 

i<k<n ( a iyfc,a2yfc,a 3 y fc ,a 4 yfc) K + ^j^^ {VkVl/tq,ty k yi/q) 

(q\ -\ y k /yi,q X k- x ty k / yi ,q K k -*i y k /yi,q K i~ K i ' y k /yi)(y k yi/q,y k yi/t) K - +K - 



i<H<n {yk/yi,yk/yi,yk/yi,yk/yi){tykyi/q,ykyi) K + +K + 

n(qh +x ty kyi /q, q K k +K ty m /q, q\ ~ X t y k /tyi){y k yi/q) x - +{ty kyi / q) x+ + 
fe Z k I 

i<k,i<n (ykvi/q, VkVi/q, yk/tq x ^yi) {ykyi/t) X7+K - (ykyi) x + +K - 

— ' — k 1 I k 1 I 

j-j. {vkVi/tq, yk/q K ' yi, Vk/tq Pl yi)xz (VkVi/q, yk/q Pl yi) x + {qyk/tq K ' yuykj 'q Kl yi) K - 

Kk,Kn 



(<? A ykvi , qyk I tq^ yi,qyk/yi) x - (*<? ft ykyi , qyk/yi) x + (<m / 1 q x ~t yi,qyk/yi) K - 

k k k 

„ {yk/q x ^yi) K + {q ai yk x i,yk/txi) x -(tq a *y k x u y k /xi) A + \ 

!<!}<„ {qyk/yi) K + ^ n {ykXuyk/tq^Xi) x -{ty k x u q-^y k /xi) x + ) ' 
l<i<m 

We give some remarks of Theorem 13. 11 As the special case n = 0, we obtain the following 
summation formula: 

E ( (-l) la]+l " +l+W Mu; Vtq/a Q ) w+Hl/ - { 

0</i-<u~<u+<ii+<a \ 
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{a 1 x i ,a 2 Xi,a 3 Xi,a 4 Xi) u - (q^ + ^ x i x j ,t 2 q^ + ^ XiXjj 

i<i< m ( t( l x i/ a i: t( l x i/ a 2,tqx i /a 3 ,tqxi/a4) u + 1<;( 4j <m (x i x j ,t 2 x i x j ) 

(q^i -Vj Xi /xj,q^ xi/xj,q u i ~ v i Xi/xj,q u * ~ v i x i /x j )(tXiX j ,qx i x j ) i 

Ki<j<m ( x i/ x j: x i/ x ji x i/ x ji x i/ x j) (t 2 XiXj ,tqX{Xj) v + +y + 



n 



(tq^ XiXj,tq Ui +Uj x i x j,q'* 1 ^ Xi/tXj)(tXiXj) - , u +(t 2 XiXj) + +v +(xi/ 'q^ xj) + 



tXiXj, Xi/tq^i Xj){qXiXj} -, - (tciXiXj) +.- {Q%i/ 
{x,iXj,Xi/q v i Xj,Xi/tq a iXj)^(txiXj,Xi/q a jXj)^+(qxi/tq u J xj,Xi/q u J Xj) v - \ 
i<i,j<m {tq a ^ +l x i x j ,qx i /tq u i Xj,qxi/xj) -{t 2 q a i +1 XiXj,qXi/xj) ^(qxi/tq^ Xj,qXi/xj) - J 
= e(u-^tq/a ) TT (tgxf , t 2 gx t 2 ) Q - -r-r (tqXiXj,t 2 qXiXj) ai 

' \a\ ///IT.- /n.i tnnr. /n.n t.nnr.j/na Irt'r; I n » \ .. A A 



h t rn (tqxi/a 1 ,tqxi/a 2 ,tqxi/as,tqxi/a4) ai 1<; ^J <m (tq a ^ +1 XiXj, t 2 q a 3 +1 x i x j ) oli ' 

(3.13) 

Note that the g-Saalschiitz sum (|2.27|) implies 



e(u; \/tg/ao)| Q |-|/3|e(-u; V g/tao)|«+|-|K-| 

r=0 V ^ r 

• V^i/ao)^!-^!^^!-^-]-^ . (3.14) 
Hence, the right-hand side of (13.121) is also expressed as 

n (Qij/fc,a2j/fc,Q3yfc,Q4j/fc)/3 fc -r-r {q Pl y k yi,tg Pl ykyi) p k rr {tyk x i,q~ a 'yk/ x i)i3 k 



x < fc < n (yl'tyDPk i<k<i< n (ykyi'tvkvfipk !< fc < n (^ ai y^,W^k 

l<i<m 



0<A-<k-<k+<A+</3 \ r=0 

-I In«-.+i_i«-n .-I I 



„ (yk/ai,yk/a2,yk/a 3 ,y k /a 4 ) K - ( q ^+K ' y kyi /tq,tq x t+ x l y m / q) 



n 



(aiy k ,a2y k ,a 3 y k ,ci4y k ) + ^} (ykyi/tq,ty k yi/q) 

(<? Afc ~ A; y k /yi,q x ^ x ' r yk/yi,q Kk ~ Kl y k /yi,q K *~^yk/yi)(y k yi/q, VkVi/t) K - +K - 

k ' i 



1<fc<Kn (Ww ' W vuVk/yu yk/yi)(tykyi/q, y k yi) K + +K + 

(q X k +x ?y k yi/q, q K " + ^y kyi /q, q X k ~ x i ' y k /tyi)(y k yi/ q) X - +K + {ty k yi/q} X + +K + 



l< k ,l<n 



(ykyi/q,y k yi/q,yk/tq x i yi)(y k yi/t) X7+K -(y k yi) : 
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{vkVi/tq,yk/q Kl yi,yk/t(i Pl yi) x-{ykVi/q,yk/q Pl yi) x +{qyk/tq K i yi,y k /q Kl yi) K -{yk/q x ^ yi) ^ 

II 



i /,./. „ {q Pl ykyuqyk/tq K ^ yuqy k /yi) x -{tq Pl y k yi,qyk/yi) x+^qyk/^ yuqyk/yi) K -{qyk/yi) K ^ 

l y k Xi,yk/tXi) x - (tq a *y k Xi, y k /xi) 

k 

(y k Xi, yk/tq a *Xi) x - (ty k Xi, q~ ai yk/xi 



n{q at y k Xi,y k /tXi) x - (tq ai y k x u y k /xi) x + \ \ 
i ; Hi 7— r 6 - • ( 3 - 15 ) 
l iii. t j iii / fn a i rA _ lt.ih.nr- a a iiii t.A . i / / 



l<i<m 

3.2 Type C case 

In the previous subsection, we derived the transformation formula of type BC by a method of 
principal specializations. In this subsection, we apply the same method to derive another type 
of transformation formula. We specialize in advance the parameters of (|2.47p so that (c, d, t) = 
(q 1/2 ,~q 1/2 ,q): 

ei=±l(i£/) i<3 i eI ° l<k<n 

^T(axf,bxf) j r (gzfa^A 



E 'II jjd 1 II /.A.±i-. 



« j_ 

/ ao<\ J. J. / St ±1\ 

Jc/ C ieJ \ « / «eJ \ x « 3 I 

8i=±i(i£j) jei c \J 



e(»; a )_ e (-i^n ^ TC rot/) n^^n 



i( qy k k / a ,qy e k /b) y\ (i 2 y e kVi l ) tt (iv^y? 1 ) 
(Vk k ) k\k Wv?) Li (yl'yt 1 ) 



K{l,..,n} \ fcGii" / fc.ieif Nafc ' k&K 

e k =±l(keK) k<l l£K c 

II 7^±T) e(u, g /a)|jcc|-|L|H -^7 11 tf ±1 

fegK Wi/fc x i / L C i^ c fcGL Wfc / k&L Wjfe 2// / 

l<i<m <5 fc =±l(A;eL) l£K c \L 



(3.16) 



Note that a = y/—ab. In this setting, the internal sum of each side simplifies drastically. 
Lemma. 3.2. 



5i=±i(ieJ) ie/ c y 



e *^n ±! ^ !fia n 7#0=(»4)' K ^ < 3 - i8 > 



, 25 k\ 11 /„,<Sfc„,±l 

'k I k£L 
5 fe =±l(fegL) leK c \L 



We can prove this lemma in the same way as in Theorem 12.11 by analyzing the residues. From 
this lemma, if u = \/— T, (|3.16j) reduces to 



v TT ^"^'K') I r ^ x i x i) TT eW) 

^ 11 11 /^ x e n 11 / fl ^ u ±i\ 

ei=±l l<i<m \ X i / l<i<j<m \ * X J ' l<i<m WX « ^fe ' 

l<j<m l<fe<n 
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e k =±l l<fe<n ' l<fc<Z<n vs, « ' l<fc<n 

l<A;<n Ki<m 



(3.19) 



Since e( % /^T;a) m _ n = (^T) m -"(-l) m ((? J"l > f a ^ 2 - , J332J is equal to 
(ax?, fee?) ir iQ 2 ^^) rr (zftff 1 ) 



(3.21) 



ETT W^j >J2-i 7 TT _ _? J / TT 
11 / 2e 4 \ 11 ( T e .i r ej ) H /r7T e '7; ±1, > 

£i=±l l<i<m \ X i / l<i<i<m \ ' J ' l<i<m W « y fc 7 
l<i<m l<k<n 

(gg^m V- TT W TT Wfi TT 9 m 

( a 2-2 m / ab ) Z_> 11 /..2e fc x 11 /yf fc yf'\ 11 / m A T ±l\ " V J ^ U J 

l<k<n l<i<m 

Let I-,I + ,K~,K + be 

I~ = {i|l < i < m, = —1}, / + = {i|l < i < m,ej = 1}, 
if- = {fc|l < A; < n,e k = -1}, if + = {k\l <k<n,e k = 1}. 

Using this notation and replacing the parameter q with q 1 ^ 2 , we obtain 

/ yy (aXi,bXi) yy -(xi/a,Xi/b) yr (qXjXj) yy (XjXj/q) 

UA J/ «> Ji- 

={1, ■■■,»*} i<j i<j 

y-r (qXj/Xj) yy _(\/M^>_\ 

jel~ l<k<n 

_ (ab) m yy \ - [ ir (y%t/g, y/qVk/b) y-r -{ayk/y/q, byk/y/q) 

~ (q l - m /ab) n \\ (Jqyk/xi) ^ 1 ^-1 ( y 2) 11 / 2\ 

l<fc<n ={l,...,n} 

yy <gg/fe2/l> TT (Vkyi/q) yy (iVk/vi) yy {\fQVkxf 1 ) \ . . 

Vi+ ^ w A i- <™> ii+ il- wvvs);" 

fc</ fc<« Zei<"~ l<i<m 

Specializing this formula as x = p a (z;q),y = pp{w;q) (a G N M ,/3 G N , |a| = m, |/3| = n), we 
obtain 

/ y-r (zi/a^j/b)^ yy (jf^Zj/Zj) yy {q^ + ^z lZj /q) 

^ \ Al (a Z i,bzi) u . Tl (zi/zo) Al (zizJq) 

0<H<a \l<i<M x 11 l<i<j<M X l/ ■ ?/ l<«<i<M \ l 31^1 



n {ziZ j /q,z i /q a iz j ) tli y-r (q h ZiW k /^/q,^/q Zi /w k ) 



{q^ZiZj^qZi/zj)^ { Zi w k / \/q, y/qZi/qP>°w k ) ^ J 
l<k<N 

(ab) H Y\i<k<N(VQ w k/a,VQWk/b)p k UiKijKM^j)^ Ui<k<i<N( w kWi) g^+ft 

Ul<k,l<N( w kWl) h 
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tt (Vqwk/q ai zi)p k ^ ( y\ ( aw k/Vq> ^WVg)^ tt {g" k Vl w k /wi) 

\<[<N <V«W*>fc ^<\li<N^ Wk / a ^ Wk l h ^l<l}l<N < Wk ^ 
l<i<M 

-|-r {q Vk+Ul W k Wl/q) yr {w k Wj / 'q,W k I 'C^ Wl) Vk -|-r {q^WkZj/ y/q, y/qW k /Zi) Uk \ 

l<i<M 

Relabeling M and N by m and n, and replacing the variables and parameters by 

Zi-> y/qxi(i = l,...,m), w k -> y k (k = 1, . . . ,n), (a,b) ->• (^/q/a 1 ,^/a 2 ), (3.24) 

we obtain a transformation formula of type C, due to Rosengren |1Q|. Corollary 4.4] . 

Theorem. 3.3. For a G N m and /3 S N™, f/ie following identity holds: 

\ - / i r (aiXj,a 2 Xj) Mi y-r (g^'M Xj/xj) yr (q^+^XjXj) 
^— ' I -J- 1 {ax;/ ai,qxi/ a.2) u (xi/xj) (xix*) 

0<V<a \l<i<m w ' ' l<«<?'<m N ' J/ l<i<j<m N 4 J/ 

-p-r : .v i .rj.x i hfj.vy :il . -p-r (q h Xjy k , qxj/y k )^ \ 

i<i}< m (q aj+lx i x j><i x i/ x j)ni (aJiy*,9 1_A ^/yfc> w / 

l<fc<n 

(g/glQ2)j«j IIl<fc<n< Q iyfc» 2yfc)/3fc ril<i,j<m<g a; i a; j)Qi ]ll<fc<;<n(^^)&+A 

(aia 2 /gl Q l)|^| rii<i<m(^i/«i> qxi/a 2 ) ai Ui<i<j<m(9XiXj} ai+aj Ui<k,i<n(ykVi) A 

tt (yk/q ai x i )p k f tt (yk/ai,yk/a2)v k tt {q Vk ~ Vl yk/yi) 

l<j<m 

tt {q Uk+Ul ykyi/q) tt {ykyi/g,y k /g^yi}u k tt {q ai yk x i,y k / x i)y k \ , g 2g . 

i<i<Kn ii,L W'vkyi'Wk/viU J k { n (yk x i,y k /q aix i}»J' 

l<i<m 

Rosengren derived this result from Gustafson's summation formula of multilateral basic hyper- 
geometric series for type C. We remark that Lassalle has derived a special case of Theorem 13.31 
from a rational function identity by the method of principal specialization [10\ Theorem 11]. His 
rational function identity (Theorem 6) corresponds to (|3.20p with a = q,b = —q. 

4 New family of g-difference operators 

In the A type case, it is known that there exists an explicit operator Hjt{u;q,t) satisfying the 
following equation [IS] : 



(«; q)ooU x A (u; q, t)V A (x; y) = (t m q n u; q)ocV y A (u; t, q)^ A (x; y), (4.1) 

where ^a(x; y) = rii<i<m T\i<k<n( Xi ~ Vk) ls the kernel function of dual Cauchy type and T> v A (u; q, t) 
is the Macdonald (/-difference operator: 

n 

V\{u;q,t)=Y J (-u) r D\M,t), (4-2) 

r=0 
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Kc{l,...,n} fceK y tf fceA' 

|Jf|=r ^ 

The operator T-C A (u;q,t) is defined by 

oo 

H x A (u;q,t)=J2u l H%i, (4.4) 



«h - E II II II T - 

l<i<j<m * J l<i,j<m yH %l 3 , H ' ^ \<i<m 

\)J.\=l 

We can obtain this fact as the special case of Kajihara's Euler transformation formula. It is also 
known that the Macdonald polynomials P A ,\(x\q,t) for type A are the joint eigenfunctions of 
H x A (u;q,t): 



Ki<m 



n A (u; q ,t)P AX ( X \q,t)=P AX ( X \q,t) [J \ utm - lq V q) ■ ^ 



The commutativity of this family {H Al }^ is proved in [21J through the Wronski relations in the 
elliptic setting. 

In this section, we give the BC type analogue of (I4.ip . Namely, we construct an explicit operator 
H x (u;q,t) which satisfies 

n x (u;q,t)^(x;y) = const. ■ V lJ (u)^(x;y), (4.7) 

n 

V y (u) = J2(-iy e (u; a) ff ,„_ r 5». (4.8) 

r=0 

Here ^(x; y) := []i<j< m e(xj; yt) is the kernel function of dual Cauchy type for type BC introduced 

l<fc<n 

by Mimachi [15] and ~ means the operation of replacing the parameters (q,t) with (t,q). Therefore, 
a = y 1 qjta and D% are the i-difference operators D% (a, b, c, d\t, q). Note that ty(x;y) is expanded 
by the Koornwinder polynomials P\(x) as follows |15| : 

*(x;y)= (-l) X *Px(x)P x *(y), (4-9) 

Ac(ra m ) 

A* = (m — A;, m - A^!, . . . , m - Ai). (4.10) 



4.1 AfRne Hecke algebras 

In order to guarantee the existence of T~L x (u;q,t) satisfying ()4.7p . we use the framework of affine 
Hecke algebras, due to Cherednik pQ and Macdonald [12]. In this subsection, we recall Noumi's 
representations of affine Hecke algebras of type C and the fundamental facts of g-Dunkl operators. 
Our notation is due to [XT] , 

We use the parameters to,t m , uo,u m which are defined such that 

II I _l ill i I _l 
(a,b,c,d) = (tZiufn, -t?nUm 2 ,q2t^u§,-q2t§u 2 ). (4.11) 
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We define the Lustig operators T§, Tf , . . . , T^j as 



t S = 4 + ^ a-^v^Ki + ^v^) ^ _ 1}j (4 . 12) 

1 — gx x 

7f = t5+H 1 ~ tei/g<+1 ( 8 f-l) (i = l > ...,m-l) > (4.13) 
II I _I 

rpx _ A i j~§ (j; ~ tmUmX m ){l + tmUm 2 X m ) (x ^ f/1 1 /H 

J m — Cm T tm " _ 2 ^ m J-J- I !. ! ! ) 

1 



Here s§, sf , . . . , s m are the simple reflections 

(*S/)0»0 = /(^r 1 ' 3T2, ■ ■ ■ , »m), (4-15) 

(47)0*0 = f{xi,---,x i+1 ,Xi,...,x m ) (i = 1,. . . ,m- 1), (4.16) 

= /(^l, *2, • • • , (4-17) 

Note that these operators sf, . . . , s m generate the Weyl group W m of type BC m . The algebra 
H(W^) generated by T§ , Tf , . . . , is isomorphic to the afhne Hecke algebra of type C m : 



(T?-t?)(T? + t i *) = (i = 0,l,...,m), (4.18) 

7?3? +1 2f = T X +1 T X T? +1 (i = l,...,m-2), (4.19) 

T?T X +1 T X T X +1 = T X +{T X T X + {T X (i = 0,m - 1), (4.20) 

If = T/lf (|i - j| > 2). (4.21) 

Here we wrote ii = ■ ■ ■ = i m _i = i. The g-Dunkl operators Yf , . . . , are defined by 

If = (If • • • T m _i)(T m • • • TJXCZ?)- 1 • • • (Iti)- 1 ) (t = 1, . . . , m). (4.22) 

We denote by C(x)[T^] the ring of (/-difference operators with rational coefficients. For any 
A x G H(W m ff ), A x is expressed as J2 w eW ^w w (^w e C(x)[T^!-]). Then we define the g-difference 
operator L X A by = ^2 w€ w A^. It is known that the following fact holds. 

For any VF m -invariant Laurent polynomial /(£) in the variables £ = . . . , £ m ), and for any 
W TO -invariant Laurent polynomial <p{x) G C[x ±1 ] Wm , one has 

/(Y>(x) = L* f{YX) <p(x). (4.23) 

Furthermore, the g-difference operator := L^ ri j satisfies for any partition A 

f{Y x )P x (x) = L x f P x (x) = P x (x)f( a t^q x ). (4.24) 

In particular, the (/-difference operators for the interpolation polynomials / = e r (£;a\t) of 
column type give rise to van Diejen's operators D x . From this view point, we call D x "column 
type" (/-difference operators. Since {e r (£; a|t)}^L 1 is the generator system of the ring C[£ ±1 ] Wm 
of W m -invariant Laurent polynomials in m variables, L x j is an element of C[Df, . . . , -D m ] for any 
/(£) G C[£ ±1 ] Wm . The operator % x {u; q, t) to be constructed in the next subsection is a generating 
function of "row type" (/-difference operators. 
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4.2 Construction of row type g-difference operators 

The results of this subsection are based on a discussion with M. Noumi. 
Let £ = (£1, . . . , £ m ) and 7/ = (771, ... , rj n ). We define 



(u; = $ (n; t) = f[ ( ^^ ± ; g)o ° , E(u; r]) = f[ e(u; Vk ), (4.25) 



where 7 is a complex number such that q 1 = t. Note that E(u; rj) is a generating function of 
e r (7y;a|t): 

n 

E(u; ri) = ^(-l) r e r (r/; a|t)e(«; a) t , n -r- (4.26) 

r=0 

This implies that 

n 

E(u;Yy)P x (y) = P x (y)He(u;at n - l q X >). (4.27) 

i=l 

Namely, the operator for / = E(u; rj) is the generating function T> y {u) of column type operators 
Dr- In the following, we regard H(u;£) as an element of u' 171 ^ C[i ±l ) Wm [[it]]. Namely, u m ~<H{u;£) 
is a formal power series in u with coefficients in the ring of W m -invariant Laurent polynomials in £. 

Lemma. 4.1. T/ze operators H(u;Y x ) and E(u;Y v ) satisfy the following identity as formal power 
series in u: 

H(u;Y x ) T , , e(u;a) n EfujY*) T , . £(u;Y») , . 

' y) = / ' j" j-T/ L ' tt(s; y) = . V ' { *{x; y). (4.28) 



H(u;atP™) y ia/ e(u;at m ) n E(u;aqP") y ' e(u;at m ) r 

Proof. Note first that E(u;aq pn ) = e(u;a) n . From f|4.9j) and (|4.24j) . the formula (|4.28p is equiva- 
lent to the identity on the eigenvalue: 

H(u;atP™q x ) = e(u;a) n E{u;agPH x ') 
H(u;atP™) e(u;at m ) n E(u;aqP") 1 ^ ''' 1 ' ' 

The left-hand side is equal to 



i<i<m (VtQuat m ~^ 1 q) Xi (y/qTtua-H- m ^q-^;q) Xt ' 



(4.31) 



Using the notation Cjj = at m % qi 1 (1 < i < m, 1 < j < n), we obtain 



(djOl) = TT ^(1 ~ y/q/iciju)(l - y/q/tCjjU l _) _ yr y/q/tCjj) 



On the other hand, since we compute 



E(u;aqPH x *) _ A { v ^uaq n - k t x t){ v / tJ^ua- 1 q k - n t- x l} 
E(u;aqP") £i {^/q~Jtuaq n ~ k ) {yftjqua- 1 q k ~ n ) 
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-q (^uaq n ~ k t x l- 1 )(^u- 1 aq n - k t x l- 1 ) 



/.. (v q/tuaq n ~ k ){^/q/tu~ 1 aq n ~ k ) 

= TT , (4.32) 

(ij)e(n")\A (^Wl/tCij) 

and hence the right-hand side of (|4.29p is equal to 

{^fqjion^ 1 ) q , n y\ (u^y/tqcjj) yr (u^ 1 ^fq/icjj) rr {u ±X y/tqcij) 

11 /„,±i nrr+ n ..\ ~ 11 /«±i/+^.,,\ 11 



□ 

Next we show that u(J^ap> m ) ^ ^ gsnerating function of the row type interpolation polynomials 
a|g, t) introduced by [8]: 

a| g ,*) = Yl rr fs™ e &> "Mb; tq^aU ■ ■ ■ e(U; t m - 1 g^ + -"^a) Vm . (4.34) 

Note that the Laurent polynomial hi(£;a\q,t) is Wm-invariant and satisfies the following interpo- 
lation property: For any partition \i 7$ (I), 

hiiat^q^alq^) = 0. (4.35) 

Lemma. 4.2. The following identity holds as formal power series in u: 

H (u;0 hi(C;a\q,t) 



"^'"^ . (4.36) 



H(u;atP™) ^ e(u;t m y/q/ia)i 
Proof. If t = q~ k (k = 0, 1, 2, . . .), from Lemma 5.4 in [8] one has 

km 

H(u;0= II e(«;^ (1-fc) Ci)fc = X)^^ a l 9 ' t ) e ^\/^/a)fcm-z- (4.37) 

l<i<m 1=0 

Since H(u;at pm ) with i = equals e(u; y/tq/a)km, by dividing the both sides of (|4.37p by 
e(u; \/ r tq/a)km, we obtain this lemma in the case of t = q~ k . In the formal power series of u in each 
side of ()4.36p . all the coefficients are the rational functions in t%. Hence the identity (|4.36p follows 
from its validity at infinitely many values of t = q~ k (k = 0,1,2,...). □ 

From this lemma, it follows that 

H(u;Y x ) ^ hiY^-a^t) 



Y, 1 '7-^ • (4-38) 



H(u;atP™) ^ e(u;t m s/q/ta)i 
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We now define the (/-difference operators Hf := Hf(a,b,c,d\q,t) (I = 0,1,2,...) to be Uj for 
/ = hi{^; a\q, t), so that 



HfP x (x) = P x {x)h l {atP™q x -a\q,t). 



(4.39) 



We call these operators Hf (I = 0, 1, 2, . . .) "row type" (/-difference operators. We also introduce a 
generating function H x (u) := H x (u;q, t) of Hf by 



H x (u;q,t) = J2 



m 



e(u;t m y/q/ta)i 

From Lemma 14.11 , we obtain a "kernel identity of dual Cauchy type". 



(4.40) 



Theorem. 4.3. The kernel function of dual Cauchy type intertwines the q-difference operator 
T-L x (u) with the t-difference operator T> y (u): 



U x (u)^{x;y) 



vy(u) 



e{u;at m ) n 

Theorem 14.31 gives the relationship between Hi and van Diejen's operators D r . 
Theorem. 4.4. For any integer I = 0, 1, 2, . . . , n, the following equation holds: 



(4.41) 



■l) l flf*(x;y)= £ 



0<s<l 



I n-l+l\ 



(4.42) 



Proof. Since HfP fl (x) = if \i C (n m ) and Z > n, 

1 



?{»v|>(x;y) = 1 + 



-HI + ■■■ + 



e{u]t m ^q/ta) e(u;t m y/q/ta) 

By using the (/-Saalschiitz sum (|2.27p . the right-hand side of (|4.4ip is expressed by 

e(u;t m a) n 



H*)y(x;y). (4.43) 



e(u; t m \J qjta)i 



£(-u 



0<r<« 



n — r 
l-r 



{t- m q'- l ,t m -'q n a z ) l ^Dy^{x;y). 

(4.44) 



Comparing the coefficient of 



e(«;t"V<?/ to )! 



in (I443|) with that in (1114]) for each Z, we obtain (|P2|) . 

□ 



4.3 Explicit formulas of iJ; 

In the previous subsection, we defined the row type (/-difference operators Hf by (/-Dunkl operators 
and showed the relationship between Hf and Df. However it is difficult to compute the explicit 
expressions of operators Hf by means of the g-Dunkl operators. In this subsection, by using the 
special case of Theorem 13. 1\ we give the explicit formulas of Hf . 
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For v G N m with < \v\ < I, we define 

(axi , bxi , cxi , dxi ) u +(a/xi, b/xi ,c/xi, d/xi ) v - 



v<i/+ \0<^-<(« m )-!/+ V l<i<m ^ilvi+v+^i fui+v, 

\u\<\u+\<l \ u -\=l-\v+\ 

n (q Vi+v iXjXj) y\ (q u ^ Uj Xi/x j ,q' / ^~ u j'xi/x j ,q^ ~^ Xj/xj) 
I ry . rv* . \ 11 / rp . I ry . ry . I ry . ry , I ry . \ 



1„-1\ 



(tXiXj,qxiXj) u + +u +(tx i Xj ,qx i x- ) v - +v - (q v i v i X{Xj, x i l Xj 



X<i<j<m^ XiX ^n+vf{ x i x 3)vf+v^ x i x j )u t +u-( x i x j ) vr + Vj \<i ,j< m {q" i+v i XiXj, q Ui+u i x i 1 x j x ) 

n(q~ v i Xi/xj) Vi (tXi/xj) +{q u ^ +1 Xj/xi,tq u t xj/xi) -(xiXj,q^ +1 x^ 1 xJ 1 ) Ui \ \ 
i « . (4.45) 

l<i,i<m ^ X il X i)vM X i/ X i) v +(q X 3/ X i^ Vi+lx 3/ X i) )) 

In particular, for |z/| = /, f G N m we have 

rr(j) r \ TT {axi,bxi,cxi,dxi) Vi yr (txjXj) Vi+V3 (q^-^Xj/xj) j-r {tXj/xj) Vi 

For any 1/ = (eii/[, . . . , e m z/^J G Z m (e« = ±1, v[ G N) such that ^^Li ^ — ^ we write |z/| = YliLi v \ 
and set H$\x; a, 6, c, d) = H^) (x^ 1 , . . . , x^ 1 ; a, 6, c, d). 

Theorem. 4.5. T/ie row type q-difference operators Hf (I = 0, 1, 2, . . .) are expressed explicitly as 

Hf= Hl l \x;a,b,c,d) J] T^. (4.47) 

i/GZ m l<t<m 
0<\v\<l 

Namely, the Koornwinder polynomials P\(x) are the joint eigenf unctions of Hf (I = 0, 1, 2, . . .): 

HfP x (x) = P x (x)h(atP™q X ;a\q,t). (4.48) 
Proof. We consider Theorem ED in the case of a = (M, . . . , M) G N m , p = (1, . . . , 1) G N n : 



I r {tqxj/ai , tqxi/a 2 , tqxj/a 3 , tqXj/a^M yr (tq M+1 XiXj,t 2 q M+1 XiXj) m 
iil m {tqxf,t 2 qxf) M i<A}<m (tqx i x j ,t 2 qx i x j ) M 



0<n-<v-<v+</i+<(M m ) V 

(01^,03^,03^,04^)^- {q^ + ^ x iXj ,t 2 q^ + ^ XiXj ) 

i<i< m ( tqXi/ ' a ^ t( l x il 'a>2,tqxi/ 'a 3 ,tqxi/ 'a 4 ) ^+ 1<;i 4j <m (xjXj, t 2 XiXj) 

n(g Mi Xi/xj,q^ ~^ Xi/xj,q u > Xi/xj,^ ~ v i Xi/xj}(tXiXj,qXiXj) v -^ 
'— 
(Xi I Xj , Xi/ Xj , Xi I Xj , X{ I Xj) (tqX{Xj , ftxiXj) + 1 + 
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n(tq^ XiXj,tq u > +u j XiXj,q^ ^ Xi/txj)(tXiXj) - + +(t 2 XiXj) + , + 



Ki,j<in (tXiXj,tXiXj,Xi/tq J Xj) (qXj,Xj) — . — (tqXiXj) + i 

n{xiXj,Xi/q v i Xj,Xi/tq M Xj) -(tXiXj^i/qMxj) +{qxi/tq u i Xj,Xi/q u i xj) 
l<i,j<m (tq M+1 XiXj,qXi/tq^ Xj,qXi/xj) ^- (t 2 q M+1 x i x j ,qx i /xj) ^+ (qxi/tq^ Xj,qxi/xj) 
{xi/q^Xj) v + {qx i yk,qXi/y k )^ 7 {tq^t Xi y k ,tq^ix i /y k ) M _ IJ + \ 

11 /,rr,/<r,\ I 11 



J 



i ,., ,,, (l x i/ x j)u+ 1<i<m {xiyk,Xi/y k ) -(tq»t+ 1 x l y k ,tqvt+ 1 x i /y k ) M _^ 

~ l<k~<n 

I rr~ I \ TT {aiy k ,a 2 y k ,a^y k ,a i y k ) T T (qy k yi,tqy k yi) 

= e{u-Vtq/a ) mM . n [[ _- ^ || 



0<A-<k-<k+<A+<(1") V 

„ (yk/ai,y k /a 2 ,y k /a 3 ,y k /a 4 ) K - ( q K+K y k y l /tq,tq x t+ x t y k y x jq) 

i<fc<n ( a iyk^ a 2yk,a3yk,a 4 yk) K + 1<k < l<n (ykyi/tq,tykyi/q) 

(q x "- x i y k /yi,q x ^~ x > r y k /yi,q Kk ~ K ' y k /yi,q K *~ K ' ~Vklm)(VkVil<l,VkVilt) K - 

:k<l 

n 



i<fe<Kn iVk/yuVk/yu yk/yi){ty k yi/q, ykyi) K + +K + 

{q x k +x ?y k yi/q, q K h +K ^y k yi/q, q x " ~ X ~ty k /tyi){y k yi/q) x - + (ty k yi/q) X + +K + 

k I k I 

[•-.a-./ „ {ykyi/q,ykyi/q,yk/tq x ^yi){y k yi/t) x - + -{y k yi) x + + - 

— ' — k 1 I k 1 ( 

n{ykyi/tq, yk/q K ' yi, yk/tqyi) x - {y k yi/q, y k /qyi) x + (qyk/tq K ^yi,yk/q K ^yi) K - {y k /q x ^yi) K + 
k k k k_ 

i k,i<n {qvkyu qyk/tq K tyu qyk/yi) x - (tqykyi, qyk/yi) x +{qyk/tq x l yi, qy k /yi) K - (qyk/yi) K + 



„ {q M y k Xi, y k /txi) x - (tq x t yk Xi, g A fe " M yfc/xj) 1 _ A + \ 

" I I " " — + + ~ • ( 4 - 49 ) 

i<k<n (Vk x i, yk/tq M Xi) x - (tq M+x k y kXi , q x k y k /xi) x _ x + J 

l<i<m 

Then the both sides are Wm-invariant for the variables x and P^-invariant for y. Since 
■j-j. (qxiyk, qxi/yk)^- Xiy k ,tq^ Xi/y k ) M _ fj + 
\<i<m {xiy k ,Xi/y k ) ^-{tq l+ ^ Xiy k ,tq»i +1 Xi/y k ) M _ + 

l<k<n 

tt T^ Xi (x,jy k ,Xi/y k ) -j-r T~^ f ~^' ) {tg M x i y k , tq M Xj/y k ) 
i<i<m ( x iyk> Xi/yk) 1 J| m (tq M Xiy k ,tq M Xi/y k ) 

l<k<n l<k<n 

the left-hand side of (|4.49p is expressed as the bilinear form: 

ST- a T T Tg* Xi {xiy k) Xi/y k ) t r T q }f ^ } (tq M Xjy k ,tq M Xj/y k ) 

Kfe<n Kfc<n 
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In terms of this form, Theorem 13.11 is regarded as a kind of bilinear transformation formula. 

We consider the special case t = q~ M in (I4.49p . Then the factors involving both variables x 
and y are expressed by 



n 

l<i<m 
Kk<n 



(qxiy k ,qxi/y k )^-(tq^ Xiy k ,tq^ Xi/y k ) M _^+ 
(xiVk,Xi/yk) ^ (tq 1+ ^ x i y k ,tqvt+ 1 x i /y k ) M _ fJ + 



n 

Ki<m 



T q %^(xf,y)T-} M -^h(x f ,y) 



^(xi;y) 



*(xi\y) 



n{q M ykXi,y k /tXi) x ~(tq x ky k Xi,q X k M y k /x i ) 1 _ x + 
k k_ 

i<k<n (VkXi, Vk/tq M Xi) x - (tq M+x t y k x, L , q x t y k /x { ) X _ A+ 



Ki<m 



n 

Kk<n 



T tm k ^{x-y k )T tm ^{x-y k ) 



$>(x;yk) 



$>(x;yk) 



(4.52) 



(4.53) 



We check that each side of (|4.49j) does not have a pole at the point t = q M . We have only to 
examine the following factor in the left-hand side: 



{q^-^/t){q- M /t) li -{q^t/t) i 
Mm ^/t){q 1 - v t/t)-{q^t/t) u 



(4.54) 



If fjJ = M for some i = 1, . . . , m, the denominator has a zero at t = q . But when \i i = 0, since 



g33D= n 



KKm 



(4.55) 



Also, when > , since (q ^/t) = (<? M /*) 



^M= I — 1 



(4.56) 



Ki<m 



Therefore the point i = q~ M in the left-hand side is an apparent singularity. Note that unless 
fi~ = or fj,f = M for each 1 < i < m, the corresponding term in the left-hand side of ()4.49p is 
zero. We can also check that the point t = q~ M in the right-hand side is an apparent singularity. 

From the argument above, specializing Theorem 13. II as a = (M m ),f3 = (l n ) and t = q~ M (M = 
0, 1, 2, . . .), we find that the factor involving both variables x and y in each side of (j4.49h simplifies 
to the form 



n 



T^{ Xi -y) 



n 



*{xi-y) 
i-.^t >0 v.iX?<M 



y(xi-,y) 



(4.57) 
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fc:A+=0 



respectively. Here we set v and follows: 



11 ilr/V 1 11 



fc:A^=l 



*(z;y fe ) 



tf(x;j/) 



(4.58) 









= M), 


'1 


(Afc = o, K 


= 0) 


z/j = < 




04<M, fi r 


= 0), K k = < 


-1 




= 1) 






V 


(otherwise) , 






V, 


(otherwise) . 





(4.59) 



In this way, the left-hand side can be interpreted as the action of a (/-difference operator on the 
kernel function of dual Cauchy type. The right-hand side is also expressed by the action of a 
t-difference operator, and in fact is equal to 



V) 1 e(u; Vtq/a ) mM -n'C y (u] a 1} a 2 , a 3 , a 4 |t, q)^(x; y). 



(4.60) 



We replace the parameters (01,02,03,04) with (a,b,c,d). Then the left-hand side of (|4.49j) can be 
expressed as 



fmM 



(4.61) 



for some o-difference operators Kf for which we will determine the explicit formulas later. Hence 
we have 



fmM 



*(x] y)~ l Yl e ( u > Vtq/a) mM -iKf (ar; y) = y) -1 e(u; V^/a) m M-n^M*(x; y). (4.62) 



Comparing (|4.41|) with (|4.62p . we obtain that 

y) _1 e(u; t m a)„e(u; y/tq/a) mM -ny- x {u)^(x; y) 

/mM \ 

= y)" 1 e(«; Vt^/a) mM -iKf V(x; y) 



,1=0 



n 

Kt<m 



(tqxi/a, tqxj/b, tqxj/c, tqxj /d) u 
{tqxj } t 2 qx1) M 



n 



(t<; M+1 rc j rc J ,* 2 <l M + 1 :Ci:c j )„ 



l<z<j<m 



( i tC£0C>i0Cj . t (£X j^X j 1 



M 



0<^"<i/-<^+<^+<(M m ) \ 



(-ir M+l " +l+l " _l e(»;^/Q) k+ |-|„-| 



n 



1 . , ^ {tqxi/a,tqxilb,tqxi/c,tqxi/d) u + v ^ 



n 

l<i<j'<m 



Xi) 

(qUi -Mj Xi j x . mi qH -Mj Xi /xj,q Vi ~ v i Xi/xj,q u * ~ v i x i /xj){tx i Xj,qx i Xj) l 
(xij Xj , x^/ , / , x^ / Xj ) (tqxiXj , t^XiXj) + , + 
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n(tq^ ^XiXj,tq u ' v i XiXj,q^ ^ Xi/tXj)(tXiXj) a - +v +(t 2 XiXj) + , 



i<j,i<m (tX{Xj) tXiXji x^/tq 3 Xj)(qXiXj)—,- (tqXiX j) + , — 

■ i j r*i j 

(xiXj,Xi/q u J Xj,Xi/tq M Xj) - (txiXj,Xi/q M Xj) + (qxi/tq u t Xj ,Xi/q u t Xj ) - 



i<!iL» (l<i u - l .r i .r j .<i.i'i/l(i"' •>',•. r/.r, /.r/) (/ 2 r/ w + l .r,.r / .r/.r ; /r ; ) (r/.r, //</''< .r,-. r/.r,/./-,-} 
{qxi/xj) u + V(x;y) 



11 (nn-Jn-A _l iDrlV */l ' ^ J 



l<i ,j<m 

Since 

e(u; i m a) n e(u; \/tq/a) mM _ n = e(u; \/tq/a) mM = H(u; at Pm ), (4.64) 
the left-hand side of f|4.63f) equals 

/mM \ 

^j/) -1 f ^e(«;V^/aW-ififl (4-65) 

Hence we have Hf^f(x;y) = Kf^f(x;y) for each Z = 0, 1, 2, . . . , mM. From the formula (I4.9P 
proved by Mimachi, we obtain HfP\{x) = KfP\{x) for any partition A C (n m ). Since n is the 
arbitrary non-negative integer, Hf equals Kf as a ^-difference operator. We see that the row type 
^-difference operators Hf(l < mM) correspond to the terms of the right-hand side of (|4.63p such 
that |^ + | — |z/ _ | = mM — I. 

We compute the explicit formula of Hf = Kf for I < M. As we will see below their coefficients 
are expressed as rational functions in frz. Note also that the operator Hf does not depend on the 
non-negative integer n. Since Hf are the Wm-invariant operators, it is enough to calculate the 

coefficients H$ '(x) := Hy\x; a, b, c, d) of Ui<i< m T qk G N m ,0 < \u\ < I). The coefficients 
Hu\x) have the following form: 

(0 



u<u- u-<v+<(M m ) 
\v\<\v-\<l \ v +\=mM-l+\v-\ 

Relabeling the indices of summations, we obtain the explicit formulas (|4.45p of the coefficients 
H$\x) for !/eN m such that < \u\ < I. 

Although we computed the explicit formula of g-difference operator Hf in the case of t = q~ M , 

for a fixed I this expression with t = q~ M is valid for any M = 1, 1 + 1, Thus the explicit formula 

(|4.45p is valid for any parameter t and we complete the proof of Theorem 14.51 □ 

4.4 Pieri formulas 

It is known that the Koornwinder polynomials have the duality property [4, 20J: 

PxiatP^q"; a, b, c, d\q, t) P^(atP™q x ; a, (3, 7, S\q, t) 



P A (o* Pm ! a, b, c, d\q, t) PJatP™ ; a, /3, 7, S\q, t) 



(4.67) 
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where the parameters a, /3, 7, 5 are defined by 



a = y abcd/q, f3 = ab/a,^ = ac/a,S = ad/a. (4.68) 

Van Diejen derived the Pieri formula of column type from the duality of the Koornwinder polynomi- 
als and D r . In this subsection, we present the "Pieri formula of row type" by using the q-difference 
operators Hi. For any partition fj, = (fj,\, . . . , fi m ), we define P+ by 

P+ = {Xe N m |Ai > • • • > X m > 0, X'j — fi'j 6 {±1,0} (1 < j < m)}, (4.69) 

where A' stands for the conjugate of a partition A. In other words, P+ is the set of partitions 
obtained by adding or subtracting at most one to /i in each column. By direct calculation, we 
obtain the following lemma. 

Lemma. 4.6. Let n be a partition. For any multi-index v E Z m , if [i+v = {^i + v\, . . . , /i m + ^ m ) ^ 

P+ , H ( u l) (atP™q>*;a,b,c,d) =0. 

For any partition /i , by substituting x = at Pm q p in (|4.48j) . we obtain 

u jg m P x {atP™) Px(atP^) 

0<\v\<l 

From Lemma 14.61 we can a PPly t ne duality of Koornwinder polynomials to (|4.70p to obtain 

£ Bg> ICff; a, t, c, rf) Y { T7/'' "t 7 

0<|^|<i 

. _ x . P u (at pm q x ;a, B, 7, <5|a. t) „„. 
P^a^a^^^g,*) 

Replacing at Pm q x and the parameters (a,/3,~/,5) with a; and (a,b,c,d), respectively, we obtain the 
following Pieri formula of row type. 

Theorem. 4.7. For any non-negative integer I = 0, 1, 2, . . . , we have the Pieri formula of row type: 



hi(x;a\q,t) - = V ff ff> (at^g"; a, ft 7, 5) , (4.72) 

lV |y ' ; P /t (a^) u j^ m 'P^ai*^ 
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